Behavior of hydrogen ions, atoms, and molecules in a-boron studied using density 

functional calculations 
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We examine the behaviour of hydrogen ions, atoms and molecules in a-boron using density 
functional calculations. Hydrogen behaves as a negative-U centre, with positive H ions preferring 
to sit off-center on inter-layer bonds and negative H ions sitting preferably at in-plane sites between 
three B12 icosahedra. Hydrogen atoms inside B12 icosahedral cages are unstable, drifting off-center 
and leaving the cage with only a 0.09 eV barrier. While H° is extremely mobile (diffusion barrier 
0.25 eV), H + and H~ have higher diffusion barriers of 0.9 eV. Once mobile these defects will 
combine, forming H2 in the interstitial void space, which will remain trapped in the lattice until 
high temperatures. Based on these results we discuss potential differences for hydrogen behaviour 
in /3-boron, and compare with experimental muon-implantation data. 

PACS numbers: 61.72.Cc,71.15.Mb,88.85.mh,13.35.Bv 
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Boron is a fascinating elemental material about which 
we are still making many new discoveries. Only recently a 
new form of Boron was found, showing unusual negative- 
U behaviour between B12 and B2 sub-units within the 
crystal 1 . Sitting somewhere between metals and insu- 
lators on the periodic table, the bonding in boron is 
strongly dependent on local environment and factors such 
as temperature, pressure and impurities. The simplest 
such impurity is hydrogen, yet to date no studies exist of 
the behaviour of hydrogen in bulk boron. 

To date 16 allotropes of elemental boron have 
been reported. The best known crystal phases are 
the a-rhombohedral structure (a-boron) and the j3- 
rhombohedral structure (/3-boron), as well as two tetrag- 
onal modifications which are probably stabilized by for- 
eign atoms^". The a-boron structure as the simplest of 
them contains 12 atoms in a rhombohcdral unit cell form- 
ing a slightly distorted icosahedron (Bi 2 , see Figure [2]). 
More complex, but thermodynamically stable at high 
temperatures is the /3-boron crystal with a unit cell con- 
taining 105 atoms 3 . Compared with the a-boron phase, 
/3-boron bulk is less dense and softer 4 . An even more 
stable modification of /3-boron has been proposed with 
106 atoms in the unit cell 5 . Boron based nanostruc- 
tures have also been proposed^, and boron single and 
double wall nanotubes and boron nanoribbons have been 
reported^—. 

Boron compounds have been proposed as promising 
candidates for hydrogen-storage 1 -, for example alkali 
doped boron sheets as possible hydrogen acceptors 1 - 1 -. Al- 
though such applied studies are already underway, there 
still remains much to understand concerning hydrogen 
behaviour in the crystalline boron phases, which forms 



the subject of the current article. Given that the B12- 
icosahedron is the basic element for both rhombohedral 
modifications, we have focused on the simpler a-boron 
phase. Based on these results comparable sites for the 
/3-boron phase are discussed. 



I. COMPUTATIONAL DETAILS 

We perform density functional theory calculations un- 
der the local density approximation as implemented in 
the AIMPRO codei^. The calculations were carried 
out using supercells, fitting the charge density to plane 
waves within an energy cut-off of 300 Ha. Electronic level 
occupation was obtained using a Fermi occupation func- 
tion with kT — 0.04 eV. Relativistic pseudo-potentials 
are generated using the Hartwingster-Goedecker-Hutter 
scheme 1 ^. These functions are labelled by multiple or- 
bital symbols, where each symbol represents a Gaussian 
function multiplied by polynomial functions including all 
angular momenta up to maxima p (I = 0, 1) and d (I 
— 0, 1, 2). Following this nomenclature, the basis sets 
used for each atom type were pdpp (B) and ppp (H) 
(a more detailed account of the basis functions can be 
found elsewhere^.). A Bloch sum of these functions is 
performed over the lattice vectors to satisfy the periodic 
boundary conditions of the supercell. Structures were 
geometrically optimized with a single k-point for the B12 
cluster in a large 23.81 A cubic supercell, and a 4x4x4 
k-point grid for pure a-boron. Both rhombohcdral and 
triclinic lattices were examined to find the optimal struc- 
ture. Hydrogen in boron was modelled using a 2 x 2 x 2 k- 
point grid for a triclinic cell containing 8 B12 icosahedra, 
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FIG. 1. Diffusion barrier (eV) for H to leave an isolated B12 
cluster. 

i.e. Bg6H„, n — 1,2. H + and H~ were modelled in charge 
neutral unit cells by the inclusion of a uniform countering 
background charge. Diffusion barriers were determined 
using the nudged elastic band method^. Saddle points 
were checked via additional energy double derivative cal- 
culations giving a single negative frequency. 

The majority of the energies considered below are rel- 
ative, for example comparison of different stable sites 
for H + within the lattice, or calculated migration bar- 
riers. There are a number of different ways to correct 
absolute binding energies to take into account interac- 
tion between charged defects in neighbouring unit cells, 
and the literature is far from being clear on the best ap- 
proach. While the standard reference for some time has 
been G. and Payne—, more recently other approaches 
are under developmen t 18 ' 19 . The base interaction energy 
can be estimated from the Madelung energy for an ar- 
ray of point charges with a neutralizing background as 
E = —aq 2 1 '(2eL), and a reasonable approximation to 
the charged cell correction is then given ini 9 - as « 2E/3. 
While corrections can be reasonably large for defects with 
higher charge states, as our H ions are charge q=±l and 
E varies with q 2 these corrections are small. Our re- 
peated unit cell sizes used for the H insertion calculations 
are also large giving L = fi 1 / 3 = 8.74A. Taking the (3- 
boron dielectric constant e «10 and an overestimate for 
the Madelung constant a of 10 only results in corrections 
to formation energies of sa0.04eV. This correction is much 
smaller than any of the energy differences we calculated 
between different hydrogen arrangements. 



II. HYDROGEN ADDITION TO ISOLATED B12 
ICOSAHEDRON 

As a basic element of the a-boron structure we first 
examine an isolated B12 cluster. Without hydrogen it 
forms a perfect icosahedron with B-B bondlengths of 
"B12 

= 1.68 A, in agreement with literature values (1.7 



AS -, 1.71 Adding hydrogen to the B12, the energet- 

ically most stable position is with hydrogen bound to the 
outside of Bi 2 (ds— h = 1-19 A), deforming the icosahe- 
dron giving bondlengths in the range cIb 12 — 1-62 — 1.74 
A. An earlier study found hydrogen stable in the center of 
the B12 cluster™, however our energy calculations (Fig- 
ure [l} show that this site is actually a metastable maxi- 
mum (with a slightly expanded B12 to c?Bi 2 = 1-72 — 1.73 
A). Any off-site motion moves the H to an interior site 
near a B-B-B triangle of the Bi 2 {cIb 12 = 1-69 — 1.83 
A). The energy difference between these two sites is only 
0.17 eV, while the position outside the Bi 2 cluster is 
6.05 eV more stable than the center position. Hydro- 
gen inside the B12 cluster migrating to the most stable 
configuration only has to overcome a barrier energy of 
AE B = 0.0042 eV. This small value indicates that hy- 
drogen atoms will sit covalently bonded to the outside of 
isolated B12 clusters. 



III. SINGLE HYDROGEN ATOMS IN a-BORON 

We next consider a-boron. For the following discussion 
we follow the nomenclature to describe inter-cage bond- 
ing adopted by-^ ' 23 i 24 , i.e. inter-layer bonds are referred 
to as 2-center (2c) bonds, intra-layer bonds as 3-center 
(3c) bonds (see Figure [2]). It should not be confused with 
the nomenclature developed for boron-hydrogen molec- 
ular compounds by Lipscomb et a@5r—, notably their 
chemical two- and three-center bonds refer to the num- 
ber of boron atoms involved in chemical bonding, rather 
than the topology of the boron cage lattice as in our case. 

Using a rhombohedral supercell we calculated lattice 
parameters to be a = 58.15° and a = 4.98 A and B- 
B bond lengths for the B12 icosahedra of d,B 12 = 1-72 — 
1.78 A with inter-icosahedral two-center bonds (2c) d,2 C = 
1.65 A and three-center bonds (3c) d,3 C — 1.98 A. The 
resulting optimised structure is in very good agreement 
with calculations from Vast et al. 23 and very similar to 
other case o 4 ' 28 . Bond lengths using a very similar triclinic 
supercell (a = 4.83 A, b = 4.83 A, c = 4.98 A, a = 60.12°, 
/9 = 60.95°, 7 = 60.97°), also with a 4 x 4 x 4 k-point 
mesh found exactly the same inter-icosahedral 2c- and 
3c-bonds lengths and d,B 12 = 1-72 — 1.77 A for the Bi 2 
icosahedra. The volume of the triclinic cell is 0.09 % 
smaller than the rhombohedral, and the total energy of 
the cell is 0.0106 eV lower. We therefore used the triclinic 
cell for all subsequent implantation calculations. 

A single neutral H atom was then added at differ- 
ent sites and the structures optimised again, this time 
with fixed lattice parameters. For H° four stable sites 
(A-D) were found, as shown in Figure [5^,b. In site A 
the H atom is located near a 2c-bond with bond lengths 
<fe A -B = 1-35 A and g?h a -b = 1-26 A. In site B the 
H atom sits centered over a boron triangle of 3c-bonds 
between the icosahedra (dH B -B = 1-27 A). Site C is 
near site B, but the H atom lies directly over a 3c-bond 
between two icosahedra (<iH c -B = 1.41 A for both). 
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As for isolated B12 we also found a much less stable 
metastable site within the B 12 cage (see TablcQI, position 
D (c?h d -b = 1-28 A for all three), where the hydrogen 
lies near a boron triangle pointing towards a 2c-bond. 
When placed in other sites within the B12 icosahedron 
(including the center site) the hydrogen moves to more 
stable locations outside the icosahedron without any en- 
ergy barrier. The energy barriers to move between site D 
and either site A or B (Table ITTI) are very small (0.13 eV 
and 0.09 eV respectively), and thus as for the isolated 
B 12 , H will not remain within the cage except possibly 
at liquid helium temperatures. This is in agreement with 
recent calculations for much higher hydrogen densities 
(1 per ~B 12 )2Q, although we find a different final stable 
location outside the cage. As well as the space at icosa- 
hedra centers the a-boron lattice also has an interstitial 
void space between the B12 layers (see Figure [2]), however 
once again hydrogen is expected to be metastable here 
and will displace to more stable neighbouring bond sites. 

We next determine the stable sites for H + and H~ 
(sites marked as A x , B x ,.. with x = +, — ). For H + sta- 
ble sites were found to be comparable to that of neutral 
hydrogen, with the most stable location near a 2c-bond 
(site A + , dn A+ -B = 1-27 A for both). Site B + is also sta- 
ble (c?Hb + -b = 1-46 A for all three bonds, directly cen- 
tered between the 3c bond triangle B atoms and inducing 
an enlargement of the 3c-bonds underneath to d^ c — 2.26 
A), but with a relative energy difference of +0.48 eV to 

(a) 



(b) 




FIG. 2. (a) stable sites B-D of H° (shown white) in the a- 
boron structure with three-center bonds (3c), (b) site A near 
a two-center bond (2c) shown in a sideview of the a-boron 
structure. 



TABLE I. Sites for H°, H+ and HT in a-boron and their 
calculated stability AE (eV) relative to the most stable site 
at that charge state (as shown in Fig. [2]). 



Site 


H° AE (eV) 


H+ AE (eV) 


H" AE (eV) 


A/A+/A- 


0.0 


0.0 


+0.95 


B/B+/B- 


+0.17 


+0.48 


0.0 


C/-I- 


+0.24 






D/D+/D- 


+2.35 


+2.31 


+3.39 


-/-/B* 






+0.56 


Inter. Void 


+0.60 


+1.82 


+0.73 



- : site is not stable in this configuration 



(a) (b) 




FIG. 3. The most stable sites for (a) H" (site B-) and (b) 
H (site B) compared with site B* for H~ over and under a 
3c bond triangle. 



site A+. For H + site C is not stable (see Table HJ). H~ 
behaves differently and this time the most stable site is 
£?_ (see Table HJ. As for H + the hydrogen sits between 
the icosahedra and induces a separation of the icosahedra 
with d 3c = 2.28 A (see Figure^). However H forms 
a shorter bond (g?h b --b = 1-2 A) and sits off-centered 
over the triangle space, pointing more into the interstitial 
void. Unlike H° and H + there is also a new site called B* 
near the B12 icosahedron on the other side of the triangle. 
On this side of the triangle the H lies over a boron bond 
of the Bi2-icosahedron with bond lengths du B ,-B = 1-35 
A and g?Hb* — b = 1-50 A (the difference between site B* 
and site B for H° is shown in Figure^). Site A_ is a 
metastable position switching ether to B_ or B*. Site C 
was not found to be stable for H~. As for H°, site D + 
and D_ (the off-centered site within the icosahedron) are 
energetically unfavourable for both H + and H _ , and we 
suppose hydrogen leaves the cage to sit in position A + 
or respectively £?_ . 

Thus in summary, the most stable site for H + and H° 
in a-boron is just off the 2c-bond between layers of B12 
icosahedra, while H _ prefers to sits in-plane between 
three icosahedra, slightly separating them. 

We next calculate diffusion barriers between the most 
stable sites for isolated H°, H + and H~ (see Table Hi]) . 
H° has a very low diffusion barrier of 0.25 eV (0.17 + 
0.08) and hence will be highly mobile in the inter-layer 
region of a-boron. H + also migrates between interlayer 
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TABLE II. Calculated diffusion barriers AE B (eV) for H°, 
H + and H~ in a-boron. (Path indicates the barrier which 
has to overcome passing from site x to site y: x — !> y) 



H° 


H+ 


H" 


Path 


AE B 


Path AE B 


Path 


A£ B 


D -> A 


0.13 


B+ -¥ A+ 0.44 


A„ ->• B_ 


0.03 


D -> B 


0.09 




B* -> B_ 


0.39 


B -> A 


0.08 









sites but with a diffusion barrier of 0.92 eV (0.48 + 0.44). 
H _ however migrates within the layers between 3c-bond 
sites, also with a barrier of 0.95 eV (0.56 + 0.39). 



IV. HYDROGEN PAIRS AND H 2 

Comparison of the total energies for two isolated H 
atoms in boron shows that the most stable configuration 
is not with two H° but rather one H + and one H , i.e. 
H in a-boron is a negative U center (see Table Mil) . This 
suggests that low density atomic H will adopt an over- 
all charge neutral mixture of H + and H - centers, with 
no net doping of the lattice. While neutral H would be 
highly mobile within the lattice, the thermodynamically 
preferred coexistence of H + and H~ centers should only 
just be mobile at room temperature. We note that, de- 
pending on the defect density and associated degree of 
electronic coupling between them via the host boron lat- 
tice, it may be possible to photocxcitc H + /H~ into a 
metastable H°/H° state, as has been shown recently for 
N2H centers in diamond^. 

While doping a-boron could encourage charge compen- 
sation by interstitial H and result in a concentration bias 
of either H+ or H - , in preliminary test calculations we 
found neither potassium or vanadium were stable within 
the B12 icosahedra. 

H + /H _ centers will experience a Coulombic attrac- 
tion and we might expect H + and H to migrate to- 
wards one another. Table IIIII shows that placing two H 
atoms in the same inter-icosahedral triangle (not dissim- 
ilar to H2 in silicon^) lowers the net system energy by 
0.84 eV compared to an infinitely separated H + /H~ pair. 
These neighbouring pairs can then combine to create an 
H2 molecule, releasing a further 0.79 eV, which sits in 
the void space between icosahedral layers (see Figure |4j. 
However it is highly constrained, and the system could 
lower its energy by a further 1.70 eV if H2 was able to 
migrate to the boron surface and escape. 

The a-boron interlayer space contains two types of in- 
terlayer void. The first and largest lies between 3c bond 
triangle sites in neighbouring layers, where the H2 sits 
(see Figure S]). These alternate with a smaller space be- 
tween a 3c bond triangle in one layer and the base of an 
icosahedron in the other. Any interlayer diffusion path 
must necessarily pass through both sites. Our calcula- 
tions show that H2 is unable to migrate as a molecule, 



TABLE III. Comparison of stability for different arrange- 
ments of hydrogen pairs within a-boron. HP refers to hy- 
drogen in its most stable location with this charge state {i.e. 
inter-cage sites). 



Hydrogen Arrangement AE (eV) 
H2 molecule separated from boron 

H2 inside the boron lattice +1.70 
H + /H~ "intimate pair" (same triangle) +2.49 
H + /H~ infinitely separated in the lattice +3.33 
H°/H° infinitely separated in the lattice +4.12 



(a) (b) 




FIG. 4. H2 in the hexagonal interstitial site between icosahe- 
dral layers, (a) side and (b) top view. 

instead dissociating into a H + /H _ pair which then diffuse 
through this smaller void before recombining. The NEB 
calculations give a total barrier of «2.2 eV, consisting of 
an rs1.4 eV dissociation energy followed by «0.9 eV mi- 
gration barriers for each hydrogen. This means H2 will be 
trapped and immobile in the large voids until high tem- 
peratures. We note that we did not find any other stable 
sites for H2 in the crystal, suggesting a maximum pos- 
sible H2 concentration in a-boron of B24H2 . The initial 
H2 separation into H + /H~ is barrierless, showing that 
H2 formation from H ions on neighbouring sites can also 
occur spontaneously. We note that we know of no other 
material where H2 is stable yet migrates by dissociation. 
The behaviour lies between that of diamond, where the 
H2 molecule is not stable, and silicon, where H2 is stable 
and can migrate as a molecule. 

V. DISCUSSION AND CONCLUSIONS 

Our results for the a-phase allow us to speculate 
on hydrogen behaviour in the more complex /3-boron 
phase. /3-boron has no equivalents of the a-boron 3c 
bond triangles, however bonds with similar lengths to 
the 2c- and 3c- bonds exist. Thus we would expect no 
equivalent of the B, B* and B~ sites for hydrogen, and 
the most stable sites are likely instead to be similar 
to A (2c-bonds) or C (3c-bonds). This agrees with 
results from implantation of positive muonium into 
/3-boron2I We can also suppose that H will not sit 
within the B12 icosahedra in the /3-boron phase. Unlike 
a-boron there are also some under-coordinated partial 
vacancy sites likely to trap hydrogen, and given the 



5 



larger connected void spaces we assume H2 migration 
will have a significantly lower barrier. 

Table IIIII shows that hydrogen within the a— boron 
lattice is less stable than isolated gas phase H2. Thus 
hydrogen incorporation within the lattice will only oc- 
cur through kinetic restrictions, i.e. if hydrogen be- 
comes trapped in the lattice, either during growth or 
through implanation, and is then not able to escape. 
Since /J— boron has larger void spaces we assume that 
H 2 will be more stable in the j3— boron void sites. Thus 
the presence of significant quantities of hydrogen within 
the a— boron lattice may serve to destabilise it with re- 
spect to the less dense /3— boron phase. This may be a 
contributing factor to the observed difficulty in synthe- 
sising a— boron, as well as the a —¥ ft— boron phase tran- 
sition (we note that molecular hydrogen in other crys- 
talline semiconductors such as silicon is often difficult to 
identify) . 

Thus we have now a first picture of hydrogen behaviour 
in pristine ev-boron. At low densities H atoms will form 
a reservoir of H + and ions in the inter-icosahedral 
space, resulting in a charge neutral crystal. Once these 



become mobile, either through thermal activation or via 
excitation into highly mobile metastable H°, they will 
diffuse together, creating interstitial molecular H 2 . This 
will remain trapped in interstitial void spaces until at 
very high temperatures when it will diffuse out. These 
results suggests that hydrogen doping of boron will be 
difficult and probably only possible through ion implan- 
tation rather than gas phase impregnation and thus a- 
boron is unlikely to serve as a useful hydrogen storage 
material. There may also be the possibility of H diffu- 
sion and resultant H2 formation at very low temperatures 
due to quantum tunnelling. 
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